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Abstract 

Recent observations of astrophysical jets emanating from various galactic nuclei strongly suggest that a double 
layered structure, or a spine-sheath structure, is likely to be their common feature. We propose that such a sheath 
jet structure can be formed magnetohydrodynamically within a valley of the magnetic pressures, which is formed 
between the peaks due to the poloidal and toroidal components, with the centrifugal force acting on the rotating 
sheath plasma is balanced by the hoop stress of the toroidal field. The poloidal field concentrated near the polar axis is 
maintained by a converging plasma flow toward the jet region, and the toroidal field is developed outside the jet cone 
owing to the poloidal current circulating through the jet. Under such situations, the set of magnetohydrodynamic 
(MHD) equations allows two main types of solutions, at least, in the region far from the footpoint. The first type 
solution describes the jets of marginally bound nature. This type is realized when the jet temperature decreases like 
viral one, and neither the pressure-gradient nor the MHD forces, which are both determined consistently, cannot 
completely overcome the gravity even at infinity. The second type is realized under an isothermal situation, and the 
gravity is cancelled exactly by the pressure-gradient force. Hence, the jets of this type are accelerated purely by the 
MHD force. It is suggested also that these two types correspond, respectively, to the jets from type I and II radio 
galaxies in the Fanaroff-Riley classification. 
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1. Introduction 

In his penetrating review paper on the active galactic nu- 
clei (AGN), Leahy (1999) stresses the importance of the di- 
chotomy of radio loud and quiet AGN and suggests the essen- 
tial difference between the respective central engines. In this 
context, our interest in the present paper is substantially re- 
stricted to the radio loud AGN and their power-down version, 
the low luminosity active galactic nuclei (LLAGN). 

It is well known (see, e.g., Urry & Padovani 1995; de 
Gouveia Dal Pino 2005) that there are two kinds of astro- 
physical jets associated with the type I and II radio galaxies 
in the Fanaroff-Riley (or FR) classification (Fanaroff & Riley 
1974). The radio emissions of the low-luminosity FR I ob- 
jects peak near the nuclei and often show symmetric two-sided 
jets. The collimations of the jets are not so strong and the 
radial velocities show considerable deceleration at large dis- 
tances. These objects are found in rich cluster environments. 
On the other hand, the emissions from the high-luminosity FR 
II objects have radio lobes with prominent hot spots and bright 
outer edges. The collimations of the jets are very sharp and 
strong asymmetries often exist between the jets and counter 
jets, with the latter being undetectable in many cases. The ra- 
dial velocities do not seem to be decelerated until the end points 
are reached. These objects exist in more isolated environments 
than FR Ts. 

Another feature we have to keep in mind is the double- 
layered structure, or spine-sheath structure, of the jets, which 
is now believed to be a common feature to the FR I and II 
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jets (Giovannini et al. 1999; Giovannini et al. 2001). The 
limb-brightened structure observed in the jets of both types 
can be well interpreted by the presence of velocity shear as- 
sociated with the double-layered structure (Giroletti et al. 
2004), in which a spine region of highly relativistic speed is 
enclosed by a mildly relativistic or non-relativistic sheath jet. 
Owing to the relativistic beaming effect, the radiation from 
the spine becomes highly concentrated within a small angle 
around the jet direction, and hence becomes weak when seen 
from other directions. The possibility that the plasma in the 
spine region may be the electron-positron pair plasma has long 
been suggested in the literature (see, e.g., Begelman et al. 
1987; Pefletier 2004). 

Historically, the problem of jet formation has been discussed 
viably in the framework of the theory of force-free or ideal 
magnetohydrodynamic (MHD) magnetospheres (e.g., Sauty et 
al. 2002; Pelletier 2004) originally developed in relation to the 
discussions of the pulsar magnetospheres. Although such treat- 
ments may be mathematically sophisticated and hence beau- 
tiful, it seems unfortunately too formal to be able to quickly 
solve many practical problems associated with the actual as- 
trophysical jets. Therefore, we adopt here a much more practi- 
cal approach to this problem. The essence of the jet formation 
and acceleration is in their non-force-free nature and non-ideal 
MHD aspects of the electrodynamic cuiTent systems. 

As stressed in our previous works (e.g., Kaburaki 2000; 
Kaburaki 2007), the jets can be formed by pinching and accel- 
erating effects of the Ampere force acting on the current car- 
rying plasma in the polar regions. The current flowing through 
the polar jet forms a returning portion of the poloidally circu- 
lating current system (see, figure 1). This current is originally 
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driven by the electromotive force caused in the accretion disc 
by the rotational motion of the disc plasma in a global, poloidal 
magnetic field. The mechanism of jet formation is thus directly 
coupled to the accretion process taking place in the same activ- 
ity site. Although there are several streams of investigations 
that attribute the formation of jets or outflows to the action of 
magnetized accretion discs (e.g., Blandford & Payne 1982), 
no satisfactory theory from the view point of the current clo- 
sure seems to exist so far in the literature. The importance of 
viewing the jet as a part of a globally circulating current system 
has been emphasized from early times by Benford (1978), and 
remarkable stability natures of such current-carrying systems 
have also been discussed repeatedly by the same author (e.g., 
Benford 2006). 

Recent observations have also revealed that both types of 
jets in FR classification are likely to be enclosed within the 
X-ray cavities, as far as they are young. The cavities are prob- 
ably maintained by the magnetic pressure of the toroidal field 
generated by the global current system circulating the disc-jet 
compounds (Kaburaki 2007). It is interesting to note that the 
shapes of such cavities may be different reflecting the differ- 
ent types of the current systems enclosed therein (see, figure 1 
and also the discussion in the final section). It is theoretically 
probable that if a jet is being driven by the accretion disc the 
resulting cavity would take an hour-glass shape, whereas the 
cavity would become of spindle shape if the accretion has been 
switched off. Observationally, the cavities of FR I jets seem 
to be of hour-glass shape with narrow necks near the equator 
(e.g., for Cen A, see figure 2 and 3 of Kraft et al. 2002, and 
figure 4 of Kraft et al. 2007; for M87, figure 6 of Forman et al. 
2007 ). On the other hand, the examples of the spindle shape 
are fairly rare. The only one plausible example hitherto known 
to the author is Cyg A (see figure 1 of Wilson et al. 2006), 
which belongs to FR II type. 

Our aim in the present paper is to construct an analytic model 
of conical MHD jets with small opening angles, which can 
reproduce the main properties of the jets described above, at 
least, in the regions well beyond their launching site. The 
method of analysis adopted here is analogous to that in our ear- 
lier work (Kaburaki & Itoh 1987) addressed to the jets in star 
forming regions. The essential difference between the situa- 
tions considered here and there is in the assumed external mag- 
netic fields. Instead of a stellar dipole field in the former paper, 
we assume here a uniform field perpendicular to the equatorial 
plasma disc, in the context of galactic nuclei. The effects of 
relativity are not included for simplicity, since the jet velocities 
in the MHD sheaths are expected to be at most mildly rela- 
tivistic. Preliminary results have been reported in Kaburaki 
(2004), though it contains the solution of only one out of two 
types obtained in the present paper 

In section 2 we first describe the global configuration consid- 
ered in the present paper, and then specify the angular depen- 
dences of relevant physical quantities in spherical polar coor- 
dinates. With these expressions, we discuss in section 3 Ohm's 
law and the equation of motion, which are essential to deter- 
mine the jet velocities. In order to obtain the solutions to this 
set of equations, we proceed in section 4 from the static force 
balance to a plausible case of dynamic balance. In the follow- 
ing two sections, two types of solutions are derived and their 




Fig. 1. A theoretically suggested classification scheme of the X-ray jet 
cavities. The paths of the globally circulating electric current are shown 
by black arrows for the case of an upward directed seed field. Also, the 
directions of the radial component of equatorial flows are indicated by 
thick white arrows, (a) Hour-glass shaped jet cavity, typical of accre- 
tion-driven jets. The appearance of a narrow neck near the equatorial 
plane is caused by the presence of inward flow towai'd the accretion 
disc. The electric current is being driven by the accretion disc around 
the central black hole. This shape seems to be a common feature of 
FR I jets, (b) Spindle shaped jet cavity, associated with a switched-off 
current system. The equatorial plasma disc is not driving the cuiTent 
system because of the lack of matter supply, and instead, the current is 
maintained by the huge self-inductance of the system. The appearance 
of the inductive electric field causes an outward flow in the equatorial 
disc (see the discussion in section 7). 

properties are discussed in detail. Finally in section 7, the re- 
sults are summarized and then a possible interpretation of these 
solutions is discussed in the light of the FR dichotomy of the 
galactic radio jets. 

2. Separation of Variables 

When the processes taking place in the central accretion 
disks are explicitly taken into account, the mechanisms of jet 
formation are expected to be as follows. In terms of the cir- 
cuit theory, the accretion disc acts as a DC current generator 
and drives poloidally circulating current system (for more de- 
tails, see section 2 of Kaburaki 2007 and figures therein). This 
is the physical cause of the toroidal magnetic component that 
is added as a twist to the originally vertical seed field. In the 
case of an upward polarity of the seed field, the radial current is 
driven outward in the disc, and dominant portion of this current 
is expected to close after circulating remote regions along the 
thin boundary layer separating the jet cavity (cocoon or hour 
glass) from surroundings. This means that the twisted field 
lines are confined exclusively within the cavity. After passing 
through the tops of the cavity wall, the currents return to the 
inner edge of the accretion disc along the polar axes. 

Under such circumstances, a polar jet can be formed from 
plasma in the polar regions, because the Ampere force exerted 
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by the toroidal magnetic field on the return current always has 
the necessary components both for radial acceleration and for 
collimation, as shown in figure 1 of Kaburaki (2007). This fact 
implies that there may always be a bipolar jet associated with 
such an accretion disc in a large-scale magnetic field, even if 
it might be strong or rather weak depending on the circum- 
stances. 

In an axisymmetric situation, the only relevant variables in 
spherical polar coordinates (r, 6, Lp) are the radial distance from 
the center r and the polar angle 6. It is often convenient to 
introduce their normalized versions, z = r/rj and rj = 9/Q, 
where rj is the footpoint radius and 8 is the half-opening an- 
gle, respectively, of a conical jet. The angular variable 77 in- 
troduced here should not be confused with that used in our 
previous papers dealing with the problems of disc accretion 
(Kaburaki 2000; Kaburaki 2001; Kaburaki 2007; in these pa- 
pers, 77 stands for (6' — 7r/2)/A, with A being the half-opening 
angle of an accretion disc). The magnetic field is generally 
written as Bq + b, where Bq is the seed field and b is the defor- 
mation due to the action of plasma motion. However, the latter 
is sufficiently large (i.e., |b| 3> |Bo|) in the accretion discs and 
jets. The influence of the seed field is taken into account only 
at the outer edge of an accretion disc (see Kaburaki 2000). 

Jets are assumed to be narrow (i.e., 8 <C 1), and 8 is re- 
garded as a smallness parameter in ordering the physical quan- 
tities. Our main interest is not in obtaining a strict solution, 
but in demonstrating the existence of a physically reasonable 
configuration that satisfies the set of equations, at least in an 
approximate sense. We therefore expect that the solution to the 
set of resistive-MHD equations can be written in a variable- 
separated form. The concrete form of each physical quantity is 
determined here by taking its symmetry property into account, 
but we do not insist that such a choice is the unique possibility. 

We adopt the following functional forms for the three com- 
ponents of the magnetic field: 



br{z,ri) = br{z) sech^r], 
bg{z,ri) = -Qbg{z)f{T]), 
b^{z,r]) = -6^(z)tanh77. 



(1) 
(2) 
(3) 



Here we have introduced the functions carrying tildes over 
them to represent the radius-dependent parts of the correspond- 
ing quantities. Although they are written as functions of nor- 
malized radius z, but they can also be regarded as functions 
of r (i.e., radius itself) whenever it is necessary. Hereafter, we 
discuss only one side of a bipolar-jet system, and hence the 
values of 6 are restricted to the range < 6' < 7r/2 whereas the 
azimuthal angle varies in the whole range < if < 2tt. 

The 7/-dependence of br is assumed to represent a localized 
structure within the cone of 6* ^ 8, which is expected to occur 
in this component owing to the balance between the sweeping 
up effect of the accreting flow at the footpoint of a jet and the 
outward diffusion caused by the presence of a finite electrical 
resistivity. The appearance of the smallness parameter 8 in bg 
indicates that this quantity is of order 0(8), and the concrete 
form of the function f{ij) is specified below from the require- 
ment of flux conservation. The presence of tanhr; in b^p guar- 
antees that bip should vanish on the polar axis and also remain 
finite outside the jet cone. 



The sign convention on the right-hand sides of the above 
expressions is adopted for an upward seed field whose r- 
component on the polar axis is positive. When the seed field is 
directed downward, Bq in the resulting solutions should be re- 
garded as having minus sign, and any other inclined case is not 
considered here. Although the seed field is amplified strongly 
by the inward motion in an accretion disc, its polarity remains 
the same. The rotational motion in the disc causes a large neg- 
ative (/3-component, which is assumed to be absent in the seed 
field. Further, the negative sign in the 6'-component means that 
a converging magnetic field in the jet structure is taken as stan- 
dard. 

For the expressions (1) and (2), the law of magnetic flux 
conservation V • b = yields 



^^i^^br) sech^T^--^^-^ 
ar rsmtl at) 



(/sin0)=O. 



(4) 



As far as we remain within the jet cone (i.e., 6 < 8), the sepa- 
ration of variable is attained approximately by the choice 



fiv) = tanh?] ^— In (cosh 77) , 

sinO 



(5) 



where the second term on the right-hand side is of the first order 
in 8 since ln(cosh 77) / sin remains finite, even for very small 
9. Then the conservation equation reduces to 



— = r— m(r br). 
br ar 



(6) 



In most cases, however, we need only the expression of f{if) 
that has the accuracy to the leading (i.e., zero-th) order in 8 as 
far as we remain in the jet-cone region stated above. In such 
cases, we can replace (5) by its approximate version. 



7(77) ~ tanh77. 



(7) 



This procedure corresponds to regarding sin in the 77- 
derivative in equation (4) as a slowly varying quantity (i.e., 
a constant) compared with a rapidly varying function tanh 77. 
Then, we can reproduce equation (6) within this approxima- 
tion. Exceptional cases in which the accurate expression (5) is 
needed are discussed in the next section. 

Given the above expressions for the magnetic field com- 
ponents, we can derive the forms of the current density from 
Ampere's law 

j = ^Vxb, (8) 
47r 

where c is the light velocity in a vacuum. Their components 
are 



d 

7^ (sine b^) 



47rrsin6' 09 
—Q^^jr{z) sech^77 

c d 



47rr dr 



9 ( , dbr 
ATrrldr^"^^' 89 



8 ^jiflyZ) sech^T] tanh 77, 



(9) 
(10) 

(11) 



where 
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jr{z) 



3e{z) 



47rrj 



Z 

z dz 

br{z) 



27rrj 



(12) 



(13) 



(14) 



The factor that expresses the order of magnitude of 
the relevant quantity has been made explicit in equations (9) 
and (11). In performing the 77-derivative in equation (9), sin 6' 
has been regarded also as a constant as noted above. In de- 
riving the last expression of equation (11), the first term in 
the curly brackets has been omitted compared with the sec- 
ond term, because the former is of order ~ 0{<d) while the 
latter is ~ O(0~^) owing to the presence of the derivative 
d/d0 = @~^d/drj. The negative sign in front of j,. means that 
the incoming current is a standard when the seed field is di- 
rected upward. It is easy to confirm that the above expressions 
actually satisfy V • j = 0. 

The functional forms of the velocity components are intro- 
duced as follows; 



Vr{z,Tj) 

V0{z,ri) = 
v^{z,ri) : 



Vr{z) sech^jy, 

-Qvg[z)Hll) 

v^{z)iBi\h.q. 



(15) 
(16) 
(17) 



Again, the hyperbolic functions and h{if) represent the exis- 
tence of a localized structure, but the sign convention here is 
somewhat different from the case of magnetic field. The jet 
proceeds outward, of course, and is expected to rotate in the 
same direction as the accretion disc and to converge toward the 
polar axis. When Vr{z), vg{z) and w<p(z) are regarded as quan- 
tities of order unity, only vg{z,r]) is of order 0(B) owing 
to the presence of a factor O. Although v^p remains finite even 
outside the jet cone, the plasma density vanishes there (see be- 
low). 

Analogously to the case of the flux conservation, the equa- 
tion of mass continuity V ■ (pv) = can be reduced to the form 

'^=r^Hr'pir), (18) 
Vr dr 

by fixing as 

Hr])^fir]). (19) 

Of course, however, some difference from the former case 
arises according to the appearance of an extra function p in 
the continuity equation. The angular dependence of the matter 
density is assumed to be 

p{z^rj)— /5(z) sech^rytanh^ r/, (20) 

reflecting a hollow-cone structure suggested by observations. 
Then, the density maximum is attained where dp/d-q — 0, i.e. 
where sech^Ty = tanh'^ yy, and we identify this surface of max- 
imum density as the middle surface of a hollow-cone jet (i.e., 
9 8). In the above derivation of equation (18), therefore, 
the term that is proportional to d(sech^ry tanh^ ?y) /drj has been 
omitted. 

The remaining thermodynamic quantities are assumed to 
take the following forms: 



T{z,i^)^T{z), 

p{z,r]) = p{z) sech^7ytanh^7y. 



(21) 
(22) 



Here, the variation of the temperature T across a jet is ignored 
since the jet is very thin. The ?y-dependence of the pressure 
represents also a hollow-cone structure. We should empha- 
size again that all of the functional forms introduced above are, 
by no means, insisted as a unique possibility for each relevant 
quantity, but that we are searching only for a plausible exam- 
ple. 

Then, the ideal gas law p = KpT is separated completely to 
yield 



p = KpT, 



(23) 



where K ^ R/p, with R and p being the gas constant and the 
mean molecular weight of the jet plasma, respectively. 

3. Ohm's Law and Equation of Motion 

The remaining fundamental equations are Ohm's law and 
equation of motion. Although the electric field is usually elim- 
inated from the scheme of MHD, especially in its ideal ver- 
sion, we treat it explicitly here. This is because it terns out 
that Ohm's law, rather than the equation of motion, plays an 
essential role in determining the radial velocity of a jet, at least 
in the region far from the footpoint. It turns out unfortunately 
that some components of Ohm's law and the equation of mo- 
tion (EOM) cannot be variable-separated completely. In such 
cases, we approximately replace the remaining angluar depen- 
dences by their representative values evaluated at the middle 
surface of a jet cone, B^^Q, where most of the jet matter is con- 
centrated. For example, whenever the dependences on scch^ry 
and tanh^ 77 appear, we evaluate their representative values as 
(sech^Ty) = (tanh^ry) = 1/2, because these satisfy the condition 
for the maximum density, sech^ry = tanh^ 77, and the identity 
sech^jy + tanh^ ?y = 1, simultaneously. 

When the system of current carrying jet satisfies a strict 
stationarity condition, the electric field should be irrotational 
and hence can be written as the gradient of a scalar potential 
U{z,rj), i.e., E = —V U. Reflecting the localized structure 
associated with the jet, we assume for the potential as 



U{z,rj) = OU{z) sech^Ty, 



and obtain 

Ee{z,T]) = 
where 



dU 

— — — = —OEr(z) sech^Ty, 
or 



im 

r~d9 
1 dU 



rsin^ dip 



Eeiz) sech^?7 tanh?y, 
0, 



Er{z) 
Eg{z) 



1 

rj dz 

2 U{z) 
rj z 



(24) 

(25) 
(26) 
(27) 

(28) 
(29) 
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The signs of Ej- and Eg here have been chosen in accordance 
with those of jr and jg, respectively. 

If the situation we are considering requires even a slow 
change of the magnetic field in time t, however, the elec- 
tric field becomes non-irrotational according to Faraday's law, 
V X E = — ( 1 /c) 9b/9i . Assuming that even in such a case the 
functional forms of each component of E given above is not al- 
tered, we can show that the only non- vanishing value remains 
in its (^-component, 

(V X E)^ = i (rEg) - 2Er\ sech^T^tanhyy. (30) 



(31) 



r dr 



Then, the (^-component of Faraday's law reduces to 



'dt 



1 d 
2d? 



rE, 



after the remaining 77-dependence has been replaced by its rep- 
resentative value at the surface of the middle plane, 9 ^ Q. 
This equation can be used to judge whether the electric field 
obtained from Ohm's law satisfies a strict stationarity condi- 
tion or not. Even if the electric field is non-irrotational, the rate 
of change for 5^ is expected to be very small since the self- 
inductance of the whole current system is very large (Benford 
2006). This point will be discussed after the explicit solutions 
are obtained. 

Ohm's law in its vector form is 



E + -V X b= 
c a 



(32) 



where a is the electrical conductivity, and is assumed to be a 
constant for simplicity. This law yields 



Er = - (veK + ^V^e) + ■ 
c (tH"^ 

Eg ^ - - [v^br+Vrbip), 

-B^ = = - {vrbg - Vgbr] 



(33) 
(34) 
(35) 



c ^92' 
when the remaining 77-dependences are replaced by their repre- 
sentative values. In deriving the above expressions, a has been 
regarded as a quantity of 0(0^^). This is because the exis- 
tence of the term jtp/a- in the (^-component equation is essen- 
tial for taking the non-ideal MHD aspects into account. As a 
consequence, the term jg /a has been omitted as small quantity 
of 0(0^) in equation (34). The origin of the resistivity may be 
classical or turbulent, but we do not need to specify it here be- 
cause a appears in the final results only through the magnetic 
Reynolds number, which we may treat as a parameter 

The MHD equation of motion in a stationary state is written 

as 

(v.V)v + i V39+-^ [bx (Vxb)]-g = 0, (36) 

where g is the gravitational acceleration, and the contribution 
from the seed field Bq has been neglected in the magnetic force 
term compared with that from the deformed part b. 

We first discuss the 6'-component of EMO. After writing out 
the full terms in this component, we keep only the leading order 
terms in the smaUness parameter 0, i.e., the terms of 0(8~^): 



K 



\ cote* : 

47r/9 I 



ld_ 



p 



Stt 



(37) 



This equation states that if we expect the magnetic pressure 
confinement of the sheath plasma (i.e., the vanishing of the 
right-hand side), then the centrifugal force caused by the ro- 
tational motion of a jet plasma should be balanced by the hoop 
stress of the toroidal field lines. Although we do not insist that 
this is the only possibility to obtain a consistent set of solu- 
tions, we examine in this paper only the cases in which this 
condition is satisfied. It is easy to suppose that such a config- 
uration would be realized in the magnetic fields developed in 
the galactic nuclear regions. Indeed, the radial component b,. is 
expected to be packed into the spine region if a converging flow 
toward the polar axis is maintained, while the toroidal compo- 
nent btp should remain weak near the polar axis as far as the 
poloidal current density is non-singular on the polar axis. 
Then, we obtain two separated equations: 



btp)^p, (38) 
where p is a function of r only. The former equation yields 



-■■2JA 



(39) 



after substitution of the representative value for 9, and the latter 
can be confirmed to hold if the condition 



P- 



K 

Stt 



(40) 



b_l 

Stt 

is satisfied. 

Similarly, retaining only the leading order terms in 9 in the 
r-component of EOM, we have 



dr 



Vg dVr 

'~~d9 
1 r, d 



r p or 



Anpr 



KQ-^i'rb^)-be — 



dbr 
~d6 



(41) 



In the above expression, the term v'^/r can be eliminated with 
the aid of the former one of equations (38). Although the grav- 
ity is usually neglected in the discussions of jet structures at 
large distances from the central object, here it is kept in the 
equation as a term of 0(1). Indeed, it turns out in the next sec- 
tion that the gravity plays an essential role in determining the 
structures of the actual jets, since the most fundamental thing 
that every jet should do is to overcome the gravity and to reach 
almost infinity. 

It should be noted that the accurate version of the function 
f{rj) in equation (5) is conceived in the definition of the direc- 
tional derivative (v • V). Then, after substituting the functional 
forms of relevant quantities and evaluating the remaining 77- 
dependences at 6* 9, we obtain 



4r V dr 

1 dp GM 
p dr r^ 



27: pr 



dlnb^f 
dr 




0, (42) 
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where A and D are newly introduced representative values, i.e., 
A = {f (rj) t'Anhr]) and D = {f{rj)/ tanh?7). We can specify the 
range of these values as0<^<l/2 and < 13 < 1, since 
/(?/) < tanhr]. However, as it is hard to specify these values 
more exactly, we treat them as parameters in this paper 

In the 95-component of EOM, all the terms survive since they 
are all quantities of 0(1). After some manipulations, we can 
reach the form 

^ _9 \ / _9 \ (43) 

^ dr r 86 J Airp \^ dr r 86 J ' 

where we have introduced the definitions 

/ = rv^ sin 6 = I sin 6 tanh 77, I = ru^ , (44) 

m = rbipSinO ~ —msinOtanhri, m = rbip. (45) 

Within the approximation that d{sm 6 tanh if)/drj ~ 
sinfl sech^r;, we have 

. -/ d\nl V0 

Vri r — ; D — 

\ dr Vr 



^ L^_Z,^Uo, (46) 
T^P \ dr br 



again after substituting representative values for the remaining 
angular dependences. 

4. From Static to Dynamic Equilibrium 

In considering the problem of jet formation, we have to care- 
fully chose the starting point from which our chain of reasoning 
to begin. As such a point, we select the static balance of a gas 
sphere around the center of gravity, apart from the transversal 
equilibrium attained in a conical jet. When neglected the mag- 
netic and inertial forces, equation (42) reduces to 

1 dp GM 
p dr r^ 



0, 



(47) 



and further with the aid of equation (23), to 
0. 



dr 



(48) 



If the temperature varies with the radius according to a 
power law. 



T(z) = Toz-" (a>0), 
then we can readily obtain the solution 

p{z) 



(49) 



PO GXp 

Poexp 



1 



-{1-. 
a I 



1} 



{a > 1) 
(a = l), 
(l>a>0) 



(50) 



where k = GAI / KT^r;]. It can be seen from this solution that 
the static equilibrium is realized (i.e., p ^ as z ^ 00) only 
for a > 1, and that the appearance of outflowing plasma is sug- 
gested for 1 > a > since then p remains finite even at z 00. 

Extending this consideration to the case of a dynamical force 
balance realized in a jet under the presence of other forces than 
the gravity and pressure gradient, we assume that p(z) takes a 
form like 



1 



(51) 



i.e., a mixture of marginal and outflow types in equation (50). 
Here, a and (3 are both constants whose values are restricted in 
the ranges, 1 > a > and /3 > 0. For this choice, we can show 
easily that 



Idp GM _ ^KT 
p dr r^ r 



(52) 



This relation indicates that, when the parameter (3 does not van- 
ish, a part of the pressure gradient force survives after cancel- 
lation with the gravity. Speaking in other way, the functional 
form of the gravity remains in equation (42) through the form 
of the miss-cancelled pressure force. On the other hand, when 
P exactly vanishes, nothing that can trace the gravity is left in 
equation (42), and only the MHD forces control the jets. 

Since the functional form of the pressure has been specified, 
it follows immediately from the equation of state (23) that 



Pqz" ^exp 



where po = po /KTq, and further from equation (40) that 

br{z) = b^{z) 



(53) 



: boz 



-1312 



exp 



(54) 



2(1 -a) 

where 69 = (87rp»o)^^^- Substituting the above expressions into 
equation (39), we have 

z}^(z) = (8i^ro)i/2^-"/2. 



(55) 



The radial velocity of the jet is essentially controlled by the 
iy9-component of Ohm's law. Substituting the expression for 
from equation (14), we can rewrite equation(35) in the from 



Vr 



27rcre2 



Woe 
Vr 



0. 



(56) 



Further inserting the results from the flux and mass conserva- 
tions, we finally obtain the expression 



-In 

dz 



br J Vr 



(57) 



where 5Rj = 2iiaQ'^r :)Voo / (? is the magnetic Reynolds number 
characterizing the jet structures. We have to note here that the 
assumption of infinite conductivity ( cr — *■ 00, i.e., the limit of 
ideal MHD) restricts the result obtained above too strictly (see 
also equation (62) below). 

Among the asymptotic solutions of the above equation (57), 
only the one that tends to a constant velocity Voo at large dis- 
tances (i.e., at z ^ 1) is of our interest here. Putting -C,. = Voo on 
the right hand side of equation (57), we obtain by integration 

(/5/Po)(Wr/uo) 



br/bo 



cxp [5Rji(l-z-i)]. 



(58) 



Further substituting expressions (53) and (54), we have 



Vriz)=VoZ^/^-" 



X exp 



2(1 -a) 



{l-z-^i-")}+5R7i(l-z-i) 



.(59) 
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In order to guarantee that this Vr (z) actually tends to a constant 
value Voo when z ^ oo, its power law dependence on z should 
vanish. Here, we have to identify this condition very carefully. 
Since the power law dependence in an arbitrary function F{z) 
causes a z^^ term in the derivative d\nF/dz, we discuss here 



dz 



InVr ■ 



Within the allowed values for a (0 < a < 1), there are two 
possibilities as for whether the second term on the right-hand 
side belongs to the first or third term: 1) when a = 1 (or 
more exactly, a ^ 1) the vanishing of z^^ term is attained 
if /3 = 2 — K, and 2) when a ~ Q the condition yields /3 = 0. 
These cases may be called the "virial-jet" and "isothermal-jet" 
solutions, respectively, because of the behaviors in tempera- 
ture. The detailed discussions of these solutions are given in 
the following two sections. 

5. Virial-Jet Solution 

5.7. asymptotic behavior 

Although we have set a = 1 in the virial jet case, actually we 
need a more careful treatment of the argument in the exponen- 
tial functions. Since a should be smaller than unity as far as an 
outflow exists at infinity, we have to write a = 1 — e {e > 0) and 
take the limit of e ^ 0. Then, the resulting expressions for the 
radial velocity is 



Vr{z) = VqZ' 



-K./2 



exp 



-k/2 



exp 



-(1 
2e ^ 

K 



1 



1 



where 



(61) 



(62) 



Since Voo represent the terminal velocity of a jet, the jet ac- 
celeration becomes the more effective the smaller the magnetic 
Reynolds number (or the conductivity) is. Especially, in the 
limit of infinite conductivity (i.e., ?Rj oo) no acceleration 
can be obtained. 

The apparent power law dependence, z^^/'^, is actually can- 
celled by the exponential factor containing e, in the limit of 
e 0. This can be confirmed with the aid of the identity 



lim 

e-+0 



z °exp 



= 1, 



(63) 



which holds for an arbitrary real number a. The constancy 
of the left-hand side of (63) is guaranteed by the vanishing of 
its derivative with respect to z, and its value is fixed as 1 by 
evaluating it at 2; = 1 . 

There is further subtlety to be considered, in this type of 
outflows. In the actual situations, every jet cannot extend to 
infinity but terminates at a large but finite radius owing to the 
interaction with the surrounding galactic or inter-galactic mat- 
ter It is therefore obvious that the jet solution too has the end 
point Tout, beyond which it cannot be applied. In such non- 
ideal cases, e can remain finite (i.e., the out flow at infinity 
does not exist) even if it may be very small. When this is 
the case, we cannot replace any exponential factors that con- 
tain e by a power law expression, by using the identity (63). 



However, as far as the jet is well developed, we may consider 
that Zout = I'ont/'r,] S> 1 and the radial velocity already ap- 
proaches a terminal velocity v^o at z = Zout- In this case, we 
have 



exp 



'2e ™* 



■^,1 ^out 



1, 



-exp^-j=l, 



(64) 



from equation (61). The former relation means that the limit of 
-Zout 00 should be taken while keeping e finite, and the latter 
equation yields 



(In Zout) ^ 



(65) 



In order to discuss the behaviors in the asymptotic regions 
(z ;2> 1) of the present kind of jets, it is convenient to intro- 
ducing a new non-dimensional radius ( = r/vout = z/zo^t in 
which the top of a jet corresponds to = 1. Then we obtain 
from equation (61) 

-k/2 



Vr{z) 



Voo C 



Z 

^out 
/2 



exp 



(66) 



where Vr{Q represents the asymptotic behavior of Vr{z) in the 
limit of z ^ Zout- This implies that the actual jets (with only 
a finite length) of virial type experience rather deceleration in 
their asymptotic regions toward their terminal values Vr{l) = 
Woo - Such a behavior is strongly reminiscent of the jets in FR I 
radio galaxies (see section 1 and references therein). 

For other quantities of our interest, the original z- 
dependences and their asymptotic forms are summarized as 
follows: 



f{z)^Toz-' 
p{z) 



:Toori = f(c), 



Po z '^^ exp 



Poo C 



PiO, 



Po z ""^exp 



-7(1-- 



hr{z) = b^iz) = boz-('-^/^^exp[-^^{l-z-^) 

— > boo r^^-'^/'^^ 6.(0 = ^(0, 

i)^(z) = (8xro)i/2;^-i/2 

= i8KTooy/^C'/'^d^{0, 



JO z-(2-«/2)exp 



joo C 



-(2-K/2) _ 



1-z- 

2e ^ 

Mo, 



2jo z-(2-«/2)exp 



-(2-K/2) ^ 



(67) 



(68) 



(69) 



(70) 



(71) 



(72) 



(73) 
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where Too = To/^out, Poo = Po/Zout' Po 



Po/zo 



bo/zont, jo = c6o/47rrj, and joo = jo/^^out- Since it is natural to 
expect that all the above quantities should decrease with radius 
C, we obtain the restriction < k < 1. 

In the asymptotic region of an actual jet of the virial type 
(i.e., e ^ 0), the continuity equations for mass and magnetic 
flux yield 



(C)=^oo( 



and 



6e(C) = 600 1 + 77 C 



— k/2 



(74) 



(75) 



respectively, since d\nVr/dC = —{k/2)(~-^ and d\nbr/dC, = 
— (1 — k/2)(^^. The (^-component of the current density is 
obtained, from the relation (13), as 

J.(C) = f JooC^'-'^/'^ (76) 

Within the jets of virial type, both sound and Alfen velocities, 
Cs and Va, vary as 



cl 



KToz-^=KT^C 



up Tip 



5.2. consistency 



(77) 
(78) 



Before the above obtained set of physical quantities are ac- 
cepted for an actual solution to the full set of MHD equations, 
some more equations are remaining to be discussed. They are 
the r- and (/3-components of EOM, and the stationarity condi- 
tion expressed in terms of the electric field. These equations 
are checked below only in the asymptotic region. 

We begin with the r-component of EOM. In the asymptotic 
region, the r-dependence of the miss-cancelled pressure gradi- 
ent force and the MHD-force are the same, i.e., —(3{KT/r) oc 
and b'^/2Trpr oc r"^ . This means that the MHD-force can- 
not dominate over the gravity even in the asymptotic region, 
and is consistent with the anticipation based on the idealized 
situation (see the previous section) that the out flow tends to the 
marginally bound case in the limit of a ^ 1 . It turns out that the 
inertial terms dominate over other forces since vf/rcc r^'^^+'^) 
and 1 < 1 + K < 2. Thus, the r-component of EMO reduces to 
the vanishing of these terms. 



c 



din Vr 



2 A 







(79) 



within the accuracy to the leading order in r-dependences, 
which yields 



A = 



(80) 



2(2 + 

For a K in the range < k < 1, we have {) < A<l/Q which sat- 
isfies the general requirement < ^ < 1/2 discussed in section 
3. 

Next, we check the (yS-component (46) of EOM. In the 
asymptotic region, / and rh behave like I ^ r^l"^ and fa ~ 



-k/2 



, respectively. Therefore, the inertial forces behave like 



dridl/dC) - Vdl/C ~ r-(i-«'/2^ and the MHD forces, like 
{br / n jj) {dm / dC) ^ bgrh/TTpC ^ r~^. In the leading order, this 
equation also requires the vanishing of the inertial terms, which 
reduces to 

D=^-^. (81) 

K + 2 

For a K in the range < k < 1, we have l/3<£'<l/2 which 
is again consistent with the general requirement < D < 1. 
Thus, in this solution, the values of the parameters, A and D, 
have been fixed consistently in terms of k. 

The poloidal components of the electric field are calculated 
from the corresponding components of Ohm's law. The 6- 
component is calculated from equation (34) as 



EeiO 



1- 



'-c\ 



(82) 



where the second term in the parentheses has been neglected 
since it is proportional to r"*^^"**)/^. From equations (33) and 
(12) we obtain for the r-component 

^ ^ ~ r j ve bg v^o A 
Er{Q = -Vrb^<— + -^ ^ + - — —z S 

C \ Vr Vr br 47rCTH^Woo Vr 



^ (1 + ^) foool'oo C\ 



(83) 



where the second and third terms in the curly brackets have 
been neglected because they are proportional to r^'^^"^/^ and 



respectively, while the first term remains finite. 



,-(1-k/2)^ 

Substituting the above obtained results for Er and Eg into 
equation (31), we have 



db^ 
dt 



2 ^ rout 



C" 



(84) 



which means that the strength of the toroidal magnetic field 
decreases with time, at least in the asymptotic region. This 
fact seems rather natural, however, because 600 — ^o/^out de- 
creases gradually with increasing Zout as the top of a jet pro- 
ceeds into the ambient plasma. Even if this is the case, the 
characteristic time for this change is estimated to be very large: 
imag ^ \d\TLb^/dt\-'^ ~ Tont/voo, which givcs tniag ~ lO"^ yr 
for rout ^ 100 kpc and Voo ^ 10^ cm s~^. Therefore, unless we 
are interested in the evolutional aspects of the jets, such a small 
change is neghgible as usually done in the current jet theories. 

5.3. global aspects 

In relation to the virial-jet solution, it should be empha- 
sized that the solution can be reconciled naturally with RIAF 
(i.e., radiatively-inefficient accretion flow)-type accretion disks 
in the external magnetic field (Kaburaki 2000; Kaburaki 
2001; Kaburaki 2007). This type of accretion model has a 
virial-type high temperature and too tenuous plasma density 
to be able to radiate efficiently. In contrast to the familiar 
RIAF models (e.g., Narayan & Yi 1994; Abramowicz et al. 
1995; Narayan & Yi 1995), however, this model explicitly 
takes into account the presence of an ordered magnetic field in 
the galactic nuclear regions. Therefore, the extraction of angu- 
lar momentum from the accretion flow is governed by the mag- 
netic stress caused by the twisted magnetic field, instead of the 
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viscous stress. It has also been clarified that a transport of heat 
energy (i.e., a non-adiabaticity) through the non-radiating disc 
causes a wind from the disc surface (Kaburaki 2001; Maruta 
& Kaburaki 2003), but this wind should be distinguished from 
jets that emanate from the inner edges of such accretion discs. 

Another feature to be mentioned is that the disc has a geo- 
metrically thin structure in spite of its high temperature. This 
is because the disc plasma is vertically compressed toward the 
equatorial plane by the magnetic pressure due to the toroidal 
field developed outside (i.e., above and below) the disc. The 
addition of this component is due to the twisting-up of the 
seed field by a rotational motion of the accreting plasma, and 
the degree of which is determined by the balance between the 
twisting motion and the slip-out (or magnetic diffusion) of the 
poloidal field under the presence of a finite electrical resistiv- 
ity. In terms of the circuit theory, the origin of this twist can 
be attributed to the existence of a poloidally circulating elec- 
tric current. As already mentioned in section 2, such a current 
is actually driven by the accretion disc operating as a current 
generator and should close finally through the polar jet regions. 

The iy9-component of the magnetic filed developed above and 
blow an RIAF-type accretion disc is given (Kaburaki 2001) by 



-{l-r 



where ^ = r/rA is the radius normalized by the size of disc's 
outer edge, ta. Further, 5Rd is the magnetic Reynolds number 
characterizing the accretion disc, and the parameter n specifies 
the strength of winds, which is a measure of the strength of the 
heat transport (or the non-adiabaticity) in the non-radiating disc 
plasma. Comparing the above expression with the asymptotic 
form {() in (70), we realize that these expressions coincide 
strictly under the identification of 

(86) 

except for the difference in the normalizing factors for the ra- 
dius. Therefore, the inequality < k < 1 reduces to < n < 
1/2, which means that there is no downward wind {n < 0) as 
far as an accretion disc is driving a jet. Substituting the expres- 
sion for K, we can specify the location of the footpoint of a jet 
in terms of the wind parameter as 



GM 



GM 



rj = > . (87) 

27iA'To KTo 

Since ?-j ^ cx) as n 0, a nearly adiabatic accretion disc 
{n ~ 0) seems to be able to drive only a weak, and therefore 
undetectable, jet. 

The coincidence of the toroidal magnetic fields generated by 
an accretion disc and by a jet implies that the origin of this field 
is one and the same current circulating through both disc and 
jet. Equating the field strength of both solutions at r = ta (i.e., 
5 = 1, or C = Ca = r A/rout < 1), we obtain 

bo = 3?D^outCi"" 1^0 1 - "^DZACr \Bol (88) 

which depends on Tout unless n = 0. This specifies the cen- 
tral field strength bo in terms of the external field |i3o|- From 
this result, we can confirm that boo = &o/^out = 5^dCa " l-^ol 
decreases with increasing 7'out since 1 — n > 0. 



For one side of a pair of jets, mass ejection rate and inte- 
grated kinetic energy flux in the asymptotic region are calcu- 
lated, respectively, as 

tt/2 



A/j(C)=/ {pvr) 2TTr^ sm9d0 
Jo 



(89) 



(90) 



Although the range of integration on 9 is extended to its full 
range (0 < 6* < 7r/2), for simplicity, this may be allowed since 
the matter density falls off rapidly to zero outside the jet cone. 
These quantities are both increasing functions of C,, reflecting 
that the velocity field has a converging component toward the 
polar axis (i.e., vq < 0). Of course, however, these values re- 
main finite with their terminal values at the top of a jet being 
Ayj(l) andi^(l). 

Three components of the Poynting flux in the asymptotic re- 



(85) gion are 



PriCv) PriC) sech^?/tanh^77, 
PeiCv) = -QPeiO sech^r^tanhyy, 



-(2-K/2) 



(91) 



PoiO = {^ErK-'^{l + l)C^'-^'\ (92) 
P^iC^v) = P^iO sech'*77tanh77, 

PAO^-^Egbr^^^C^'-''/'''. (93) 
47r 47r 

The radial flow of electromagnetic energy through one side of 
a jet pair is defined as 



S{0= I 2TTr'^PrSmede = -e^rjblvoo c 

6 



k/2 



(94) 



This quantity is again an increasing function of (, reflecting the 
convergence of the magnetic field (i.e., bg < 0). 
The ratio of the Poynting flux to the kinetic flux. 



(95) 



is a decreasing function of and its terminal value is 
S'(l)/F(l) = (80/3)(A'To/u^Zout)- If we assume roughly 
that the dominance of Poynting flux should be kept in the virial 
jets all the way to their end points, the terminal velocity can be 
estimated as 



80 KT( 







^oo- , , (96) 

•J ^out 

by setting S'(l)/F(l) ~ 1. In this case, the terminal velocity 

— 1/2 

decreases with the jet length according to z^^^ , which also 
reflects the marginally trapped nature of the jets of this type. 
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6. Isothermal- Jet Solution 

6.1. consistency 

The isothermal-jet solution is specified by the set of param- 
eters, a = and /3 = 0. In this case, the relevant quantities take 
the following forms: 



f (z) To const, 

p[z) =poexp{-K(l-z"^ 

p{z) = poexp{-K(l - 



)} =]5ooexp(KZ ^) 

)} =Pooexp(Kz~^) 



br{z) 



he{z) 

^oexp|-|( 



1-, 



Vr{z) = woexp 



■)} 



{( 



fid-.- 



v^{z) = (Si^To)'/-' = const., 

where poc, = Po exp(-K), poo = Po exp(-K), 
6oexp(— k/2) and — tjq cxp(k/2 + 
Since 

-ln4„ = -l„V = -2= . 



^1 ~ 

— mVr = 

dz 

— Ini 
we obtain 

be{z) = 6o 



z'^pVr) = 



)^{2-(f-*-' 



= 6„ 



2- -z" 
2 



cxp < 



CXp( -2 



-')} 



and 



X exp ^' 
= i;oo{2-(| 



3?- 



fid-. 



2 



X exp 



Further, for the characteristic velocities, we have 



p _ 


KTa 


p 






= % 


np 


TTp 



8A'Tn 



: const. 



(97) 
(98) 
(99) 

, (100) 

(101) 

(102) 

600 = 



(103) 



(104) 



(105) 



(106) 

(107) 
(108) 



In contrast to the virial-jet solution, the isothermal-jet solution 
has an essentially exponential structure, and the quantities re- 
main finite even at large distances (z 00). 



As in the case of the viral-jet solution, we have next to 
check the consistency of the isothermal-jet solution. In the 
r-component of EOM, the inertial term tends to Av'^ since 
Vr Voo as z cx), while the MHD-force term tends to 
—SKTqD. Therefore, in the leading order, the equation re- 
duces to 



(109) 



rv^ cx r and 



vi,A - 8KT0D = 0. 

As for the (y3-component of EOM, we find I 
rh = rbip oc r in the asymptotic region (z ^ 1). Since the 
inertial forces and MHD forces become of the same orders 
of magnitude, i.e., {br/Trp){dm/dz) ~ bgrh/npz ~ z° and 



Vr{dl/dz) 
dlnl 



D 



vel/z' 



, equation (46) is held only when 



dlnm 



dz Vr dz 

and both of them are satisfied if 

Further, if the approximation 

fiv) (/(^)tanh7y) 



br 



2A 



(110) 



(111) 



(112) 



\tanh77/ (tanh^Ty) 

is allowed, we can fix as ^ = 1/4 and D = 1/2, and the terminal 
velocity of a jet of this type is specified as 



z;oo=4(/vro)i/2. 



(113) 



Speaking more generally, it can be concluded from equations 
(109) and (108) that the flow becomes super- Alfvenic as far as 
D> A. 

The poloidal components of the electric field are obtained 
from Ohm's law. In the asymptotic region (z ^ 1), they are 



brVr 



1 




(114) 



and 



Er{z) 




(115) 



where the third term in the curly brackets has been neglected. 

As in the case of the viral jets, these field components do not 
satisfy the stationarity condition. The associated change rate of 
the magnetic field is given by 



db^ 
'dt 




(116) 



The characteristic time for the changing magnetic field is given 
by ^mag ~ roat/voo, which is also large as in the case of the 
virial jets. In both cases, the long time scales can be understood 
as a consequence of large inductances of the relevant current 
systems (Benford 2006). In the equations above and in the the 
next subsection, we may substitute the value in (1 1 1) for D. 
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6.2. global aspects 

Since we do not know any analytic solution of the accre- 
tion disk in a large-scale magnetic field, which is naturally 
connected to the isothermal-jet solution, we cannot specify the 
strength of the magnetic field bo, in relation to the accretion 
process. The mass and kinetic energy flux through one side of 
a jet are calculated in the asymptotic region as 



An 



M,{z) = — Q'rip^ 
F{z) = YQ^'r]p^vl 



(117) 
(118) 



These are increasing functions of z, since the jet flow is again a 
converging one. Although these quantities diverge as z oo, 
practically there is no jet that extends to infinity. If we want to 
know the exact nature of a jet expanding into real vacuum, then 
we must trace its entire evolution. However, such a problem is 
beyond the scope of the present paper 

The three components of the Poynting flux in the asymptotic 



region are 



47r 




Pe{z) 



P^{z)^^^Egbr 

in 



Then, the integrated flux becomes 



Sriz) 




(119) 



(120) 



(121) 



(122) 



This is also an increasing function of z, reflecting the conver- 
gence of the magnetic field (i.e., bg < 0). The ratio of the 
Poynting flux to the kinetic flux is given by 



S{z) 
F{z) 




(123) 



which is independent of z. Thus, the isothermal jets seem to 
remain always Poynting-flux dominated. 

7. Summary and Discussion 

We have found two types of jet solutions to the set of MHD 
equations, which are valid in the region except very close to the 
footpoint, under the restriction that the jet plasma is confined 
within a thin sheath structure centered on a conical surface of 
small opening angle. Commonly to both types of jet, such a 
density enhancement of jet plasma (assumed to be a normal 
plasma consisting of electrons and ions) is maintained in the 
valley of the magnetic pressures that is formed between the 
region of poloidal-component dominance near the polar axis 
and that of azimuthal-component dominance outside the cone. 
In this configuration, further, the centrifugal force due to jet 
rotation is always balanced by the hoop stress of the toroidal 
magnetic field. Although we have not discussed explicitly for 



it, the spine region is very likely to be filled with highly rel- 
ativistic electron-positron pair plasma as various observations 
seem to suggest. The terminal velocity of a jet of both types 
of solution does not exceed by far the thermal velocity at its 
footpoint, and so remains to be at most mildly relativistic. 

The virial-jet solution is characterized by a virial-like tem- 
perature profile, and it seems that the jets are marginally bound 
by the gravity. Namely, although the radial velocity of a jet 
increases toward a terminal velocity at infinity in the case of 
free expansion, when the jet terminates at a large but finite 
distance, as is the case in every actual situations, the velocity 
decreases according to a power law in the asymptotic region 
to reach a terminal value. All other quantities also vary ac- 
cording to power laws. This type of solution can be matched 
quite naturally with the RIAF type solution of accretion discs 
in a large-scale magnetic field (Kaburaki 2000; Kaburaki 
2001; Kaburaki 2007). This fact strongly suggests that the 
virial jets are driven in the actual situations by such accretion 
dicks. The ratio of the Poynting flux to kinetic flux decreases 
with the radius, and may decrease to a value smaller than unity. 

On the other hand, the isothermal-jet solution is character- 
ized by a temperature independent of the radius, and the jets are 
accelerated purely magnetohydrodynamically, since the grav- 
ity is exactly cancelled out by the pressure gradient force. The 
radial velocity is gradually accelerated and finally seems to 
reach a super-Alfvenic value. In this solution, however, the 
quantities other than the temperature vary as exponentials of in- 
verse radius, and the changes in their values from the footpoint 
to the top remain rather small. The ratio of Poynting to kinetic 
flux is a constant in radius and remains always Poynting-flux 
dominated. Further, there is no accretion disk solution hitherto 
known to us, which can be matched naturally with this type of 
jet solution. 

The correspondence of the virial jets to PR I jets may be 
already evident from the above description. The FR I jets are 
known from observations to have well defined opening angles 
and decreasing radial velocities at large distances. The latter 
fact is explained here by their marginally bound nature. The 
opinion that the matter entrainment is the cause of deceleration 
(as often suggested in the literature, e.g., Bicknell 1994) is 
not confirmed in our solutions. Although the virial-jet solution 
indeed has the flow and magnetic field that converge toward the 
jet axis, and hence the so-called entrainment of surrounding 
matter is realized there, such a configuration is also found in 
our isothermal-jet solution in which the radial velocity does 
not experience a deceleration. 

The host of FR I jets are ellipticals often in rich cluster en- 
vironments, and their accretion discs are believed to be radia- 
tively inefficient. The smooth matching of the virial-jet so- 
lution to the RIAF-type accretion disc solution seems to be a 
strong support of the idea that a RIAF in a global magnetic 
field always drive a virial jet. The strength of such a jet is very 
likely to be controlled by the parameter k, which specifies the 
height of jet's footpoint. Since this value has been determined 
as twice of the wind parameter n (defined in Kaburaki 2001), 
which is a measure of non-adiabaticity in the accreting matter, 
the strength of a virial jet may also be controlled by the strength 
of such a radial energy transport taking place in the accretion 
disc. 
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Although the connection of the isothermal-jet solution to FR 
II jets may be not so clear as the connection of the virial-jet 
solution to FR I jets, it seems nevertheless to be possible. The 
overall trend that the radial changes of various quantities are 
rather mild in this solution seems to be favorable to FR II jets. 
In contrast to the virial-jet solution, the jet velocity in this solu- 
tion increases gradually until the end point is reached, in spite 
of the presence of an entrainment (i.e., a convergent flow). 

Finally, we refer to a rather speculative possibility of inter- 
preting FR II jets, i.e., the idea that FR II jets are dying sys- 
tems, in which their central engines have been switched off. 
Originally, this idea was proposed with a support found in the 
trend seen in the power-size (P-D) diagram of FR II sources, 
which seems to suggest an decrease in the jet luminosities with 
their ages (Baldwin 1982). If this is the case, a jet that has 
been driven once as a RF I jet by the gas supply owing to some 
past event of galaxy encounter is now switched off, for the lack 
of matter supply, and is dying as a FR II jet. Also, we may see 
recurrences of the jet phenomena in rich clusters, because there 
may be periods of inefficient matter supply between successive 
galaxy encounters. In fact, examples of such recurrences can 
be seen in some examples as the successive bubbles seen in 
X-rays (for M87, see e.g., Forman et al. 2007). 

Other kind of radio sources, such as compact symmetric ob- 
jects (or CSO, see e.g., Readhead et al. 1996) that resem- 
ble FR II sources but are much smaller and hence believed to 
be young, have also been interpreted as dying systems. The 
problem of over population that may arise when we assume an 
evolutionary track from CSOs into FR II sources can thus be 
avoided. Most of the young sources will die young as CSOs 
(e.g., Giroletti 2000). Only a small fraction of young systems 
can evolve as FR I sources to a big size, and then die as well- 
developed FR II sources when their central engines (or more 
properly speaking, dynamos) are switched off. Therefore, large 
FR II sources such as Cyg A seem to be very rare cases. 

If we take this interpretation of FR II sources seriously, it 
becomes likely that the morphology of their jet cavities, which 
are revealed by recent X-ray observations, are different from 
those of FR I sources. Jet cavities of RF I sources seem to be 
typically of hour-glass type that have narrow necks near the 
equatorial plane. This fact can easily been understood as re- 
flecting the presence of infalHng matter associated with their 
accretion discs (see figure la). On the other hand, the shape 
of the jet cavity of a switched off accretion system is likely to 
have a spindle shape with no neck in the equatorial plane (see 
figure lb). This is understood as follows. The poloidal current 
flowing through a switched-off jet is maintained by the self- 
inductance of the system. The current flowing in the jet closes 
its circuit through the equatorial plasma disc and through the 
walls of the X-ray cavity in the northern and southern hemi- 
spheres. Different from the accretion-driven cases, however, 
the inductive electric field appears in the equatorial disc having 
the same sign as the electric current, trying to keep the current 
constant. Since the magnitude of the electric current cannot 
increase in spite of this inductive field, the rotational motion 
of the disc is forced to increase in order to satisfy Ohm's law, 
equation (33). The unbalanced centrifugal force caused by this 
super-Keplerian rotation drives the disc plasma outward, re- 
sulting in a secretion disc (an inverse of the accretion disc). 



In this respect, it is interesting to note that Cyg A, a typical 
of FR II sources, seems to have a spindle shaped X-ray cavity 
(Wilson et al. 2006) although further observations are needed 
for confirmation. The equatorial rings of hot gasses seen in 
their X-ray images might be reflecting the existence of out- 
flows. Anyway, if the above quoted interpretation of FR II jets 
is correct, similar spindle shaped X-ray cavities are expected 
also for other FR II sources. 
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